DUDLEY  WWX  LIBRARY 
(MAL  POSTGRADUATE  SCHOOI 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 

DATA 

REDUCTION  FOR  THE  UNSTEADY  AERODYNAMICS 
ON  A  CIRCULATION  CONTROL  AIRFOIL 

by 

Billy  Murel  Pickelsimer 

March  1977 

Thesis  Advisor:            Louis  V.  Schmidt 

Approved  for  public  release;  distribution  unlimited. 


178054 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whmn  Dmtm  Inttr.d) 


DUDLEY  KNOX  LIBRARY 
NAVAL  POSTGRADUATE  SCKOOt. 


REPOftT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

1.    REPORT  NUMBER 

2.  GOVT  ACCESSION  NO. 

3.     RECIPIENT'S  CATALOG  NUMBER 

4.    TITLE  (mnd  Subtltlm) 

Data  Reduction   for   the  Unsteady 
Aerodynamics   on   a  Circulation 
Control  Airfoil 

5.     TYPE  OF  REPORT  *  PERIOO  COVERED 

Master's   Thesis 
March   1977 

•  ■    PERFORMING  ORG.  REPORT  NUMBER 

7.     AUTHOR^ 

Billy  Murel   Pickelsimer 

t.    CONTRACT  OR  GRANT  NUMBERS; 

».    PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

Naval   Postgraduate   School 
Monterey,    California        93940 

10.     PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  m  WORK  UNIT  NUMBERS 

II.    CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Naval   Postgraduate   School 
Monterey,    California        93940 

12.     REPORT  DATE 

March   1977 

13.     NUMBER  OF  PAGES 

61 

14     MONITORING  AGENCY  NAME  *   ADDRESSfJ/  dlltmrmnt  from  Controlling  Ollleo) 

Naval   Postgraduate   School 
Monterey,    California        93940 

IS.     SECURITY  CLASS,  (ol  thl*  rdport) 

Unclassified 

IS«.     OECLASSIFICATION/OOWNGRAOING 

SCHEDULE 

16.     DISTRIBUTION  STATEMENT  (at  unit  Kmrnort) 

Approved   for   public   release,    distribution   unlimited. 

17.     DISTRIBUTION  STATEMENT  (ol  (ha  mbmtrmct  mntmrmd  In  Block  30.  II  dlltmrmnt  from  ftoport) 

18.     SUPPLEMENTARY  NOTES 

19.    KEY  WORDS  (Contlnuo  on  tmrmrtm  ml  mm  tl  nmcmmmmrr  mnd  Identity  by  mloek  numbmt) 

20.    ABSTRACT  (Cmnttnum  on  tmrmrmm  »ldm  II  nmcmmmmry  «*4  Identity  by  blmmk  mmubmr) 

Calculating   the    lift,    drag,    and  pitching   moment    coefficients 
for   an   airfoil    from   the   static  pressure   distribution   obtained 
from  wind  tunnel    tests    is    routine    task  when   steady   flow   is 
considered,    but    it    is   much  more   complicated  when   the    airfoil    is 
operating    in   an   unsteady   flow   field,    similar   to   that   experienced 
by   a  helicopter   rotor  blade,    produced  by   an   oscillating   wind 
tunnel.      A   data   reduction   routine   capable   of   condensing   the    large 

DO  ,;Srn  1473 

(Page    1) 


EDITION  OP  1  NOV  SS  IS  OBSOLETE 
S/N   0102-014-6601  ! 

1 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OP  THIS  PAOE  (9hmn  Dmtm  Kntmrmmt) 


UNCLASSIFIED 


fbCU*»lTV   CLASSIFICATION   OF   THIS  PtflEf^w   D»i«  Ent»r»d 


numbers  o£  data  associated  with  the  unsteady  investigation, 
as  well  as  a  numerical  integration  algorithm  for  the  unsteady- 
aerodynamic  coefficients,  were  developed;  however,  no  unsteady 
data  were  collected  due  to  hardware  failures.   The  ability 
of  the  program  was  demonstrated  on  previously  obtained  steady 
and  quasi-steady  data  and  sample  results  were  presented. 


DD  Form   1473 
„,  1  Jan  73 
S/N  0102-014-6601 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  F»AGErW»>«"  Dmtm  Enffd) 


Approved  for  public  release,  distribution  unlimited 


Data  Reduction  for  the  Unsteady  Aerodynamics 
on  a  Circulation  Control  Airfoil 


by 


Billy  Murel  Pickelsimer 
Lieutenant,  United  States  Navy 
B.S.,  Oklahoma  State  University,  1970 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  AERONAUTICAL  ENGINEERING 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
March  1977 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


Rear  Admiral  Isham  Linder  Jack  R.  Borsting 

Superintendent  Provost 


This  thesis  prepared  in  conjunction  with  research 
supported  by  the  Naval  Air  Systems  Command  370C. 

Reproduction  of  all  or  part  of  this  report  is 
authorized . 


Released  as  a 
Technical  Report  by: 


ABSTRACT 

Calculating  the  lift,  drag,  and  pitching  moment  coeffi- 
cients for  an  airfoil  from  the  static  pressure  distribution 
obtained  from  wind  tunnel  tests  is  routine  task  when  steady 
flow  is  considered,  but  it,  is  much  more  complicated  when  the 
airfoil  is  operating  in  an  unsteady  flow  field,  similar  to 
that  experienced  by  a  helicopter  rotor  blade,  produced  by 
an  oscillating  wind  tunnel.   A  data  reduction  routine  capable 
of  condensing  the  large  numbers  of  data  associated  with  the 
unsteady  investigation,  as  well  as  a  numerical  integration 
algorithm  for  the  unsteady  aerodynamic  coefficients,  were 
developed;  however,  no  unsteady  data  were  collected  due  to 
hardware  failures.   The  ability  of  the  program  was  demonstrated 
on  previously  obtained  steady  and  quasi-steady  data  and  sample 
results  were  presented. 
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I.   INTRODUCTION 

The  primary  objectives  of  this  investigation  were: 
1)  to  report  on  the  data  reduction  techniques  used  to  ob- 
tain values  for  both  steady  and  non-steady  aerodynamic 
coefficients,  resulting  from" the  experimentally  acquired 
static  pressure  distribution  over  the  model  when  the  tunnel 
was  operating  at  a  specified  dynamic  pressure,  with  Coanda 
blowing  oscillated  at  various  frequencies  and  amplitudes  of 
momentum  blowing  coefficient;  and   2)  to  present  a  sample 
of  results.   Removal  of  the  model  from  the  tunnel  for 
modifications  in  November  1976  precluded  any  actual  process- 
ing of  non-steady  data  consistent  with  the  stated  goals. 

However,  the  data  reduction  techniques  are  applicable 
to  the  pressure  distribution  obtained  with  steady  Coanda 
blowing  and  will  be  demonstrated  on  such  data  obtained  in 
April  1976.   Prior  to  removal  of  the  model,  the  tunnel  was 
operated  without  blowing  in  an  effort  to  generate  a  set  of 
non-steady  data  to  be  used  in  tailoring  the  data  reduction 
procedures  to  be  used  downstream  in  the  on-going  project. 
In  that  test,  the  model  experienced  self-excited  harmonic 
pressures  resulting  from  vortices  being  shed  alternately 


The  modifications,  to  correct  for  surface  delamina- 
tions  near  the  upper  jet  nozzle,  include  a  steel  upper 
surface  trailing  edge  and  slot  height  adjustment  to  0.015 
inches.   Additional  pressure  taps  were  installed  on  the 
upper  surface  forward  of  the  Coanda  slot,  and  leaking 
tubes  were  repaired  while  the  modifications  were  being 
implemented. 


from  the  upper  and  lower  surfaces.   The  manually  recorded 
data  consisted  of  the  RMS  and  mean  values  of  the  unsteady- 
pressure,  and  the  phase  angle  in  respect  to  a  reference  for 
each  pressure  orifice.   A  more  complete  description  of  the 
data  acquisition  system  is  given  in  Section  II  B. 
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II.   EQUIPMENT  AND  INSTRUMENTATION 

A.   AIRFOIL 

The  airfoil  used  in  this  study  was  obtained  from  the 
Lockheed-California  Company  after  completion  of  their 
design  feasibility  study  on  circulation  control  rotors 
(Ref.l).   It  had  an  elliptical  cross-section  with  a  chord 
of  10.215  inches,  thickness  ratio  of  0.20,  and  camber  of 
3.3%  (Figure  1).   A  slot  to  facilitate  Coanda  blowing  was 
machined  in  the  upper  surface  at  chord  position  X/C=0.95, 
with  the  slot  height  set  at  0.010  inches.   Coanda  air  was 
supplied  to  the  plenum  cavity  at  both  ends  of  the  section 
to  insure  uniform  flow  along  the  span.   Fifty-three 
pressure  taps  at  mid-span  with  a  chord-wise  orientation  were 
monitored  through  two  scanivalve  transducers  for  an  elec- 
trical representation  of  the  pressure  distribution  over  the 
airfoil.   The  pressure  tap  locations  are  listed  in  Table  I. 

The  2.0  x  2.0  foot  unsteady  flow  wind  tunnel  in  the 
Department  of  Aeronautics  Laboratories  at  the  Naval  Post- 
graduate School  was  used  as  it  is  capable  of  simulating  the 
velocity  flow  field  experienced  by  blade  elements  of  a 
helicopter  in  forward  flight.   The  model  was  mounted  at  an 
angle  of  attack  of  -5  degrees,  which  closely  corresponds  to 
the  zero-lift  angle  of  attack  for  the  data  reported  herein. 
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B.   INSTRUMENTATION 

The  data  acquisition  system  evolved  from  a  classic 
manometer  board  arrangement  with  readings  hand  recorded 
(Ref.  2)  and  reduced  to  the  present  system,  which  employed 
electrical  read-outs  with  an  existing  Digitec  system  with 
customized  signal  conditioning.   The  digital  multiplexing 
circuitry  in  the  Digitec  system  was  interconnected  through 
the  scanivalve  solenoid  controllers  for  automatic  scani- 
valve  advancing  as  the  output  of  the  digital  voltmeter  was 
sequentially  printed  on  a  paper  tape.   Pressure  sensing  was 
accomplished  by  connecting  each  of  the  pressure  taps  to 
one  of  the  two  scanivalves  (Figure  2)  through  tubing  of 
uniform  diameter  and  length.   This  was  done  to  insure  a 
constant  transfer  function  for  all  pressure  taps,  when 
measuring  unsteady  pressures.   Calibration  tests  using  the 
method  reported  by  Johnson  (Ref. 3)  showed  an  acceptable 
variation  in  dynamic  gain  in  the  frequency  range  from  0  to 
300  Hz  (Figure  3) ,  which  included  the  range  of  practical 
interest. 

For  the  unsteady  aerodynamics  analysis  (C  =0) ,  the 
voltage  representation  of  the  surface  static  pressures  was 
of  the  form 

ex=ax+bx(t).  (2.1) 

Full  analog  measuring  techniques  were  employed  in  measuring 
the  mean  voltage,  RMS  voltage,  and  phase  relationship  to  a 
reference  for  each  pressure  tap.   The  mean  voltage  was  ob- 
tained by  processing  the  signal  through  a  low-pass  filter 
circuit  (Figure  4)  with  a  two  second  time  constant  which 
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corresponds  to  a  corner  frequency  of  0.08  Hz  on  a  Bode  plot. 
Therefore,  the  influence  of  the  unsteady  voltage  was  essen- 
tially negated  and  the  output  readings  represented  the 
steady  components  of  the  signal. 

C.   ADVANCED  DATA  ACQUISITION  SYSTEM 

To  meet  the  primary  objective  mentioned  in  Section  I, 

a  more  sophisticated  data  acquisition  and  reduction  system 

2 
will  be  required  and  is  now  being  constructed  by  Englehardt 

at  the  Naval  Postgraduate  School  in  the  Department  of  Aero- 
nautics.  The  system,  based  on  the  Intel  MDS-80  micro- 
processor unit,  will  control  the  logic  of  the  experiment, 
collect  and  store  data,  and  handle  the  entire  data  reduction 
process.   The  algorithm  for  data  reduction  discussed  in 
Section  V  was  proven  sufficiently  simple  and  capable  of 
performing  all  required  data  reduction  and  numerical  integra- 
tion. 


2 
Master  of  Science  thesis  by  Englehardt  with  a  complete 

description  of  the  entire  microprocessor  system  will  be 

published  in  June  1977. 
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III.   DATA  REDUCTION 


A.   STEADY  DATA  PROCESSING 

The  aerodynamic  coefficients  are  defined  by  integral 
equations  but  must  be  approximated  by  numerical  integra- 
tions, since  data  are  available  only  at  a  finite  number 
of  pressure  tap  locations.   The  approximate  steady  normal 
force,  chord  force,  and  pitching  moment  coefficients  are 


1.0 


Cn  = 


0 


VVd(  yc> 
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Y/C  (max) 


C   = 
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'm(TE) 


Y/C  (min) 
Y/C  (max) 


VVd(/^} 


(3.1) 


-Y/C  (min) 
1.0 


CCp£-cpr)(Y/c)d(Y/c)  + 


0 


VV(1"X/c)d(X/c) 


Then,  with  an  angle  of  attack  correction  for  the  force 
coefficients  and  a  moment  transfer  to  the  quarter-chord 
position,  the  aerodynamic  coefficients  are 


C- -C  cosa-C  sina 
In      c 

C,=C  sina+C  cosa 
d  n      c 


Cm(c/4)"Cm(TE)"0,75Cn- 


(3.2) 


The  pressure  data  from  the  April  1976  tests  and  the 
November  1976  mean  pressure  data  were  transferred  to  a 
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Hewlett-Packard  9830  tape  cassette  and  processed  using  the 
program  listed  in  Appendix  A,  which  employs  the  equivalent 
numerical  relationship  to  the  equations  stated  above. 

B.   UNSTEADY  DATA  PROCESSING 

Before  the  unsteady  data  were  suitable  for  numerical 
integration,  a  dynamic  gain  correction  was  necessary  to 
account  for  the  transfer  function  discussed  in  Section  II. 
The  frequency  of  the  self-induced  oscillations  was  constant 
at  262  Hz  when  the  tunnel  airspeed  was  approximately  100 
ft/s.   From  Figure  3,  the  corresponding  gain  correction, 
|G-| ,  was  approximately  0.5,  which  implied  that  the  amplitude 
of  the  surface  static  pressure  was  attenuated  by  a  factor 
of  one-half  at  the  transducer.   All  phase  angle  values  were 
relative  to  tube  number  7,  which  was  arbitrarily  chosen  as 
a  clock  reference.   Phase  angle  errors  resulting  from  the 
two  different  scanivalve  circuits  were  faired  out  in  the 
analysis  (Figure  5). 

With  the  above  corrections ,  the  numerical  integration 
of  the  unsteady  static  pressure  distribution  was  carried 
out  in  a  manner  similar  to  the  steady  pressure  integrations. 
For  the  unsteady  normal  force  coefficient 

,  1.0 

/      2v  *  f  X 
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Then 


And 


u  -i  l&  J- 

<Cn2>2   =    ^|(<Cn2>2sin^2+(<Cn2>2cos*) 

*T  =   tan"1      (<Cn2>sin*> 
(<C    z>cos*) 


2   ** 
where      (C      )        is    the   RMS  value   of   the  normal    force   coefficient 
x  n   ' 

and   <j>T      is    the  phase   angle   relative   to   tube   number   7.      The 
chord   force   and  pitching  moment   coefficients   were   calculated 
similarly. 
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h 


.    (<C/>  sin*) 

T  s — 

(<Cc2>  cos*) 


/     2\  2      . 
<Cm  >     sin* 


Y/C(max) 


te 


-Y/C(min) 


(C      sin*-C      sin*)(Y/f)d(Y/    )    -f 
Pf  Pr  u  u 


1.0 


(C     sin*-C      sin*)(l-X/r)d(X/r) 
pl  pu  l  L 
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2    z 
<C      >     cos  +- 


m 


te 


Y/C(max) 


Y/C(min) 


(Cp   cos<j>-Cp   cos*)(Y/c)d(Y/c)      + 


(C     cos<f)-C      cos4»)(l-A/r)d(A/r) 
Pi  Pi- 


u 


«-t- 1 


(<Cm2>  sin^2+(<Cm2>cos*)2      ,   *T 


m 


m 


(<Cm2>2sin<D) 

— T^ — 

(<Cm   >  cos*) 


The  unsteady  aerodynamic  coefficients  are  calculated  as  in 
Eqs .  3.2,  with  the  aid  of  elementary  trigonometric  identities 
to  account  for  the  phase  angles.   The  program  listed  in 
Appendix  A  has  the  option  for  unsteady  numerical  integration 
incorporated. 

It  is  intended  that  the  C  =0   oscillating  pressure 

y 

measurements  be  repeated  when  the  improved  data  acquisition 
system  becomes  available.   The  tests  will  be  conducted  at 
several  tunnel  dynamic  pressures  to  determine  the  Reynolds 
number  dependence  of  the  data.   The  previous  measurements 
were  at  approximately  R  =0.54  x  10  . 
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IV.   RESULTS 

A.   STEADY  DATA  RESULTS 

The  validity  of  the  concept  of  circulation  control  was 
demonstrated  by  Englar  (Ref s .  4  and  5)  and  is  confirmed 
here  with  the  steady  momentum  blowing  coefficient  results. 
The  resulting  aerodynamic  coefficients  are  summarized  in 
Table  II.   The  ability  to  produce  relatively  large  in- 
creases in  lift  coefficient  with  only  a  small  increase  in 
blowing  coefficient,  shown  graphically  in  Figures  6a,  b, 
and  c,  is  one  of  the  desirable  characteristics  of  the 
system.   Probably  the  best  measure  of  the  effectiveness  is 
the  aerodynamic  lift  amplification  ratio,  3C,/3C  ,   which 
was  computed  to  be  approximately  15.3,  the  low  value  being 
attributed  to  the  geometric  deficiencies  of  the  trailing 
edge  slot  as  alluded  to  in  Section  I.   However,  the  results 
are  encouraging. 

Close  scrutiny  of  the  coefficient  of  pressure  informa- 
tion plotted  against  X/C  (Figures  7-12)  reveals  a  pressure 
distribution  over  the  model  such  that  the  center  of  pressure 
was  at  the  mid-chord  position,  which  is  in  agreement  with 
the  theory.   The  center  of  pressure  position  is  also 

verified  by  the  graph  of  C  ,  ,„>.  Vs .  C  . 
3  6  v  m(c/4)      y 


(h  =  £   -  CW   =  .25  -  ("3-5)  =  .48 

L  C.P.   4    (3C../3C  )         (T3TT) 

1    y 
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TABLE.  II 
AERODYNAMIC  COEFFICIENT  SUMMARY 
STEADY  AERODYNAMIC  COEFFICIENTS 
Angle  of  Attack=-5  degrees  U=100  ft/s 


Run  No. 

C^ 

Cl 

Cd 

Cm(c/4) 

42301 

0 

-0.0942 

0.0341 

-0.1065 

42302 

0.0855 

1.0332 

0.1083 

-0.3907 

42303 

0.0640 

0.8146 

0.0678 

-0.3365 

42304 

0 . 0448 

0.5774 

0 . 0440 

-0.2776 

42305 

0.0272 

0.3459 

0.0327 

-0.2175 

43206 

0.0133 

0.1100 

0.0367 

-0.1664 

C1=-0.0765 


cd=  0.0451 

C  =-0..1063 

m 


UNSTEADY  AERODYNAMIC  COEFFICIENTS 

<f>=225  degrees 
</)  =  39  degrees 


C1(RM3)=2.1404 


Cd( RMS) =0.1629 
C  (RM3)=0.5817 


<f)  =  63  degrees 
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Also,  the  low  drag  qualities  of  the  airfoil  are  favorable 
when  full  thrust  recovery  of  the  momentum  blowing  coefficient 
is  assumed.   Pressure  coefficient  data  are  listed  in 
Appendix  B. 

B.   UNSTEADY  DATA  RESULTS 

The  unsteady  aerodynamics  of  the  elliptically  shaped 
airfoil  have  been  given  a  rather  cavalier  treatment  in 
the  past  but  their  contribution  to  the  total  forces  and 
moments  must  be  realized.   The  implications  of  the  unsteady 
aerodynamics  is  not  generally  well  understood;  therefore, 
a  simple  example  is  presented  here  to  emphasize  their 
importance  before  discussion  of  the  unsteady  results. 

Consider  an  airfoil  at  zero  angle  of  attack  with  the 

RMS  pressure  distribution  equal  on  the  upper  and  lower 

surfaces  (Figure  13a)  but  with  unequal  relative  phase 

angles  (Figure  13b)  taken  from  the  quarter-chord  and  the 

three-quarter-chord  positions  from  the  actual  unsteady  data 

(Figure  5).   The  unsteady  normal  force  coefficient,  \C  )    , 

is  calculated  from  the  numerical  approximation  to  Eq.  3.3, 

as  : 

(c      >   sin<J>  =    0.5(sin(210°)+sin(230°)-sin(35O)-sin(550)) 

/     2^ 

(C      >  sin<t>    =    -1.33 

x  n    ' 

(C      )  cos*    =   0.5(cos(210O)+cos(230°)-cos(35°)-cos(550)) 

(C      )  cos<}>    =    -1.45 

<Cn2)"2   =     ^jl.332    I    1.452'  *    =    I^ff 

(CnZ)      =    1.97  4>    =    223° 
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Similarly,  for  the  moment  coefficient  (referenced  to  the 
trailing  edge) : 

(C  2>  sin*  =  0.5((sin(210o)-sin(35°))0.75+(sin(230o)-sin(55°))0.25) 


m 


2  2 

<C  >  sin<j>  =  -0.6 
x  m  '     Y 


2.* 


<C   >  cos  =  0.5((cos(210O)-cos(35°))0.75+(cos(230O)-cos(55°))0.25) 


*s 


<Cm2>  cos  =  -0.78 


i* 


<Cm2>   =  >.62  +  0.782      ♦-   ^ 

?    % 
(C      >  =  0.98  $  =    217° 

x  m  / 


Clearly,  the  force  and  moment  coefficients  are  of  significant 
magnitude  even  if  a  uniform  distribution  (RMS)  of  pressure 
coefficient  were  assumed,  because  of  the  importance  of  con- 
sidering relative  phase  relationships.   The  dynamic  behavior 
of  the  unsteady  aerodynamics  of  the  circulation  control 
airfoil,  particularly  the  lift  amplification  ratio,  in  the 
unsteady  environment  caused  by  cyclic  blowing  will  be  in- 
vestigated after  the  improved  data  acquisition  system  is 
operational. 

As  discussed  in  Section  I,  measurements  were  taken  in 
November  1976,  with  the  airfoil  experiencing  harmonic 
pressure  variations.   The  unsteady  data  are  listed  in 
Appendix  C.   The  experimental  data  were  numerically  inte- 
grated, similar  to  the  previous  example,  for  the  lift,  drag, 
and  pitching  moment  coefficients  (Figures  14  and  15) .   It 
should  be  recognized  that  the  phase  shift  between  the  upper 
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and  lower  surfaces  is  approximately  180  degrees,  the  result 
being  the  relatively  large  amplitude  of  the  aerodynamic 
coefficients.   Highlights  of  the  unsteady  results,  summarized 
in  Table  II,  were  that: 

1.  The  time  varying  pressures  over  the  model  were 
predominantly  harmonic,  with  a  non-dimensional  frequency 
(Strouhal  no.)  of 

S  =  t^   =  0.449. 

2.  The  RMS  value  of  the  normal  force  coefficient  was 
approximately  2.14,  which  was  an  order  of  magnitude  greater 
than  the  RMS  drag  coefficient. 

3.  The  test  Reynolds  no.  was  approximately 

R  =  0.54  x  106. 
e 

4.  The  location  on  the  airfoil  for  the  center  of 
pressure  due  to  unsteady  loadings  was  approximately  0.491 
chord  and  the  RMS  value  of  pitching  moment  at  this  reference 

point  was 

2  h 
(C  L)      =  0.1694 
x  m 

with  a  phase  angle  of  90  degrees  lag  relative  to  the  normal 

force.   This  term  may  be  considered  as  an  apparent  mass 

type  effect. 


36 


-Z.Hr 


CPCRM5)*5lN(PHl)  V5.  X/C 


UPPER  5URFHCE 


-1.0. 


+       + 


+     + 


♦   t 


+ 
* 


0.0^t 


y 

a 

III 

in 
y 

a: 
a. 


1.0. 


0.2 


0.4 


0.E 


0.8 


^  X/C 

1.0 


L 
P 


2.0L 


z 
y 

-  -w 


L 

L 

y 
a 


-LB. 


LOWER  SURFACE 


in 

z 


0.0 


+       0.2 


H.H 


0.6 


0.B 


1*0 

♦ 

* 


x/c 


LI. 


2.1 


Figure    14 


-2. EL 


-1.0. 


a. a 


1.0. 


u 
z 

3 
ill 

in 
y 

Q. 
L 

a 

i- 
z 
y 

-    -2.1 


2.0L 


CPCRM5)KD5CPHl)  V5.  X/C 


UPPER  5URFRCE 


+       +       + 


+  ♦ 


*+ 


0.2  0.H  0.G  0.B 


h= d ^->y< 


L 
k 

y 
a 

IT: 

z 


-1.0. 


0.2 


2.0L 


LQWER  SURFACE 


0.2  e.H  0.G  0.B 


Bfr        *    +    *    x  * 


1.0 


*         + 


Figure    15 


38 


V.   DATA  REDUCTION  FOR  THE  ADVANCED  DATA  REDUCTION  SYSTEM 

The  advanced  data  acquisition  system  described  in 
Section  II  C  will  allow  a  more  comprehensive  investigation 
of  the  unsteady  aerodynamics  of  the  model,  but  the  data 
reduction  techniques  must  be  tailored  to  handle  the  large 
numbers  of  data  to  be  collected.   Oscillation  of  the  plenum 
pressure  as  described  by  Bauman  (Reference  6)  induces 
harmonic  blowing  from  the  trailing  edge  slot  at  a  specified 
frequency,  W..  .   It  is  assumed  that  the  unsteady  static 
pressure  at  each  station  will  have  a  steady  value  plus  a 
periodic  value  with  frequency  W-  superimposed  due  to  the 
oscillating  momentum  blowing  coefficient. 

Since  the  data  are  assumed  periodic  in  form,  the  un- 
steady pressure  at  any  point  on  the  airfoil  may  be  defined 
by  a  Fourier  series  (Ref.  7)  as  follows: 

P(x,t)=A0(x)+n|1  [An(x)  cos(n.W1t)+Bn(x)  sin^t)]    (5.1) 

where  W..  is  the  frequency  of  the  oscillation.   The  standard 

formulae  for  obtaining  the  Fourier  coefficients  for  the 

infinite  series  are: 

T 
P(x)dx 
-T 
T 
AR(x)  =  i    P(x)cos(nx)dx  (5.2) 

-T 

f  T 

P  (x)sin(nx)dx 
-T 
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A  (x)  =  -i 

0      2T 


Bn(x)  =  J 


Because  the  experimental  data  are  in  the  form  of  a 
discritization  of  a  truncated  periodic  signal,  the  above 
equations  are  slightly  altered  for  estimating  the  Fourier 
coefficients.   The  modified  form  is: 

*       1   K 

£  (x)  =  i  z   p  (x) 

o      K  ±ml      i 

*        2  K 

A  (x)  =  ^  E   P. (x)cos(nW-t)   ;   n=l,...,n    (5.3) 
n      K  i=1  i        1 

2   K 
B  (x)  =  £     Z   P. (x)sin(nw  t) 
n      K   i=1   i         1 

where  k   is  the  number  of  samples  to  be  analyzed  and  n 
is  the  harmonic  number.   It  is  recognized  that  the  physical 
meaning  of  A  (x)  is  the  steady  or  mean  value  of  the  oscillat- 
ing pressure.   That  is,   A  (x)  =  (jp    (x))   where 


<p„00>* 


V  O 


Lim   1_ 
o^J/       T+  co   2T 


P(x)dx  (5.4) 


By  subtracting  out  the  mean  pressure  value,  the  signal  becomes 

biased  to  a  pure  oscillating  pressure  with  a  mean  value  equal 

a,  ^ 

to  zero.   The  terms  A  (x)  and  B  (x)  are  equivalent  to  the 

n         n 

cosine  and  sine  components,  respectively,  discussed  in 
Section  III  B.   Figure  16  illustrates  the  reconstruction  of 
a  discritized  signal  utilizing  this  algorithm  programmed  on 
the  Hewlett-Packard  9830  calculator  (program  listing  in 
Appendix  A) .   The  error  was  minimized  by  insuring  that  the 
number  of  samples  analyzed  was  an  integer  multiple  of  the 
frequency;  however,  the  truncation  error  was  never  greater 
than  2%. 
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Having  developed  an  effective  method  for  handling  the 
large  numbers  of  data,  calculation  of  the  aerodynamic  force 
and  moment  coefficients  can  be  accomplished.   Before  the 
actual  processing  of  data  begins,  the  sample  pressures  are 
converted  to  non-dimensional  coefficient  form  to  expedite 
calculations.   The  unsteady  pressure  coefficients  are 

Cp-(t)  =  H^H      i  =  1»2»3>---»53     C5.5D 

where  P  and  H  are  steady  pressure  values  from  a  pitotstatic 

tube  mounted  on  the  floor  of  the  wind  tunnel.   Then  the 

oscillating  signal  at  any  pressure  tap  is; 

C  (x,t)  =  I   A   cos(nw1t)+B   sinfnw^t) 
pv  *  ^      n    vl'n    K      1    J 

C  (x,t)  =  I  T   sin(nw.t  +  <j)  )  (5.6) 

p  n    v   1   Yn  ^   ' 


-i    2    2  B 

where         T  -  \A  +B    .   $   =   n, 

n   |  n   n    '   yn       /. 

An 
The  unsteady  C  (x,t)  are  related  by  the  phase  angles   $  , 
however,  the  phase  angles  are  meaningless  unless  they  are 
measured  relative  to  a  "clock"  or  reference  signal.   There- 
fore, it  is  necessary  to  designate  one  pressure  tap  to  be 
used  as  the  reference.   Then,  the  relative  phase  angle  is 

Trel   yi   yref  v   ' 

For  simplicity,  <\> .      shall  hereafter  be  interpreted  as  the 
relative  phase  angle. 

With  the  steady  pressure  coefficient  removed  from  the 
oscillating  signal,  separate  aerodynamic  coefficients  can 
be  calculated  for  the  steady  and  unsteady  signals.   The 
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steady  force  and  moment  coefficients  are  found  exactly  as 
in  Section  III  A.   The  unsteady  aerodynamic  coefficients  can 
be  defined  in  a  Fourier  series  form  as: 

Cn(t)  =  E  Cn   sin(kw1t  +  (|)k) 

Cc(t)  =  Z  Cc   sin(kw1t+({)k)  (5.8) 

C  (t)  -  E  C   sinfkw.  t+<b,  ) 
mv        m,     v   1   Tky 

The  magnitudes  and  phase  angles  of  the  coefficients  can  be 

determined  thusly: 

for  the  first  harmonic 

C      sinfw,  t  +  <j>)    =   C      sin<|>cos  (w..  t) +C      cos*sin(w..  t) 
n,         v    1      YJ  n,  v    1    '      n,        Y         v    1 


where 


C     sin*   = 

nl 
1  0 


fc      (-)sin<J..  (-)-C      (-)sin*    (£j|  d(£) 
P-,    c  Yl^c^      p    ^c  ru\CJ      VCJ 


C      cos* 

nl 


C      (£)  cos  <>-(-)    -    C      (-)cos(j)    (-)!   d(£)       (5.9) 
p,    c  Tl   c^  p      c  yuvc^       KcJ       v        J 

0        l  u  J 


then 
C 


j = ^  ,    .C      sin*, 

n      =VCT1    C0S*iJ      +    «"„   sin^r,        ♦      =   Tan"1  /    nl  L 

nl        \     nl  -1  nl  -1  -1  1  C     cos?* 


nl 


Similar  calculations  are  done  to  obtain  values  for  the  higher 

harmonics  (C   ,C   ,  ' *  *  C   )  as  well  as  the  chord  force  and 
n2  n3»      ni 

pitching  moment  coefficients  using  the  higher  harmonics  of 
pressure  coefficient  information-  as  obtained  during  the 
initial  data  reduction  previously  described. 
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VI.   CONCLUSIONS 

As  mentioned  in  Section  I,  the  two-fold  objective  of 
this  paper  was  to  develop  a  computer  program  for  unsteady 
aerodynamic  analysis  and  to  demonstrate  its  capabilities 
on  experimentally  generated  unsteady  data.   The  objective 
was  partially  satisfied;  a  computer  program,  well  within 
the  limitations  of  mini-calculators  such  as  the  Hewlett- 
Packard  9830,  was  developed  which  employed  a  straight- 
forward numerical  integration  of  pressure  coefficients  to 
obtain  values  for  the  aerodynamic  coefficients.   The  pro- 
gram can  be  modified  for  unsteady  data  by  properly  taking 
into  account  the  cosine  and  sine  components  of  the  pressure 
signal  relative  to  a  fixed  harmonic  (clock)  signal.   Also, 
the  extraction  of  the  Fourier  coefficients  of  the  truncated, 
time-varying  signal  appears  to  be  within  the  capabilities 
of  the  updated  data  acquisition  system. 
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APPENDIX  A 
COMPUTER  PROGRAM  LISTINGS  AND  OUTPUT 


1  REM  PROGRAM  FOR  DECOMPOSITION  OF  A  FUNCTION 

2  REM  INTO  FOURIER  COEFFICIENTS 
10  DIM  EC2561 

30  W1=(20*PI) 

35  DISP  "DIVISIONS  PER  CYCLE?"; 

36  INPUT  D 

42  P=2*PI/W1 

43  T=P/D 

50  REM  COMPUTE  SAMPLE  FUNCTION 

60  FOR  1=1  TO  220 

65  X=W1*(I-1)*T 

70  E[I>SIN(X)+0.5*SIN(2*X)+0.3*SIN(3*X) 

80  NEXT  I 

100  PRINT 

105  PRINT  "THE  DRIVING  FREQUENCY  IS  10  HZ" 

106  PRINT  D;  "DIVISIONS  PER  CYCLE" 

107  PRINT 

110  PRINT  "     COEFFICIENTS" 

111  PRINT 

113  PRINT  TAB2"COS"TAB15"SIN" 

114  FOR  J=l  TO  3 

115  PRINT 

117  PRINT  TAB2"A";J,  "B";J,  "PHASE  ANGLE" 

119  A1=B1=0 

120  FOR  1=1  TO  220 

130  Bl=(Bl  +  ECn*SIN(J*Wl*(I-l)*T)) 

140  Al=(Al+ECIl*COS(J*Wl*(I-l)*T)) 

150  NEXT  I 

159  P=ATN(B1/A1) 

160  WRITE  (15,170) (Al*2/220) ,  (Bl*2/220)  ,  (P*180/PI) 
170  FORMAT  F6 . 3 , 7X,F6 . 3 , 13X,F5 . 1 

175  NEXT  J 

180  END 


THE  DRIVING  FREQUENCY  IS  10  HZ 
10  DIVISIONS  PER  CYCLE 

COEFFICIENTS 

COS         SIN 

A  1         B  1  PHASE  ANGLE 

0.000       1.000  90.0 

A  2         B  2  PHASE  ANGLE 

0.000       0.500  90.0 

A3         B  3  PHASE  ANGLE 

0.000       0.300  90.0 
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1  REM  COMPUTE  CL , CD, CM- -STEADY  AND  UNSTEADY 

2  REM  FROM  CP  DATA  ON  TAPE  CASSETTE 
10  COM  R1,DC80J  ,P[801  ,C£60] 

20  DIM  AC1061 
30  PRINT 

40  LOAD  DATA  4, A 
60  C9=0 
70  X9=0.95 
80  Y9=0.048 

90  DISP  "CP  DATA  FILE  NO.  IS"; 
100  INPUT  Z 

110  LOAD  DATA   Z 

111  DC223=    0.4 

112  DC36J  =  0.539 

113  Gl=l 

115  PRINT  "INPUT  OPTION" 

116  PRINT  "1  FOR  STEADY  DATA" 

117  PRINT  "2  FOR  NON-STEADY  DATA" 

118  INPUT  G3 

119  GOTO  G3  OF  130,212 

120  GOTO  130 
130  PRINT 

140  PRINT  "RUN  NO.  IS"R1 

150  PRINT  TAB7"UPPER"TAB24"LOWER" 

160  PRINT  TAB2"N"TAB9"CP"TAB19"N"TAB26"CP" 

170  FORMAT  F3 . 0 , 2X,F7 . 3 , 5X,F3 . 0 , 2X,F7 . 3 

180  FOR  N=l  TO  27 

190  M=N+26 

200  WRITE  (15,170)N,C[N]  ,M,C[M] 

201  NEXT  N 

202  FIXED  3 

203  PRINT  "CP1="C9"AT  X/C="X9",  Y/C="Y9 

204  STANDARD 

205  DISP  "CORRECT  THE  CP  VALUES" 

206  STOP 

210  GOTO  270 

212  DISP  "INPUT  DYNAMIC  GAIN"; 

213  INPUT  G2 

214  FIXED  3 

216  PRTNT  "RUN  NO.  IS"R1 

217  C9  =  C[2l]  +  0.33*(C[22]  -C[21]) 

218  PRINT  "TAP"TAB7"CP (MEAN) "TAB17"CP (RMS) "TAB2  7"PHI" 
220  FORMAT  F3 . 0 , 3X,F7 . 3 , 3X, F7 . 3 , 3X,F4 . 0 

222  FOR  1=1  TO  53 

224   WRITE    (15, 220)1, C[IJ ,D[I J ,P [I] 

226  NEXT  I 

227  PRINT  "C9="C9 

228  DISP  "CORRECT  THE  CP  VALUES" 
230  STOP 

235  GOTO  270 

240  FOR  1=1  TO  53 

242  C[I]=G2*D[I]*SIN(P[IJ/57.296) 

244  NEXT  I 

2  45  C9=(D[21]  +  0.33*(D[2  2]  -D[21]  )) 

246  Gl=2 

247  GOTO  270 
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250    FOR   1=1   TO    53 

252    C[IJ=G2*D[I]*COS(P[I]/57.296) 
254   NEXT    I 

256  C9=(D[21J+0.33*(D[22]-D[21])) 

257  Gl=3 

258  GOTO    270 
270    L1=0 

272  L2=0 
274  P1=0 
276  P2=0 
310    FOR   1=1   TO   26 

320    IF    1=21   THEN   490 

330   Cl  =  -C[rj 

340    C2  =  -C[I  +  1J 

350    X1=A[IJ 

360    X2=ACI  +  1] 

370  GOSUB  1120 

380    L1=L1+N1 

390    P1=P1+M1 

400    Cl  =  C[I  +  27] 

410   C2=C[I  +  26] 

420    Xl=A[I  +  27] 

430    X2=A[I  +  2$| 

440    GOSUB    1120 

450    L2=L2+N1 

460    P2=P2+M1 

470    NEXT    I 

480    GOTO    640 

490    Cl  =  -C  [21] 

500    C2=-C9 

510    X1=A[21] 

520    X2=X9 

530    GOSUB    1120 

540    L1=L1+N1 

550    P1=P1+M1 

560    Cl=-C9 

570    C2=-C[22] 

580    X1=X9 

590    X2=A   22 
600    GOSUB    1120 
610    L1=L1+N1 
620    P1=P1+M1 
630    GOTO   400 
640    D1=0 
650    D2=0 
660    P3=0 
670    P4=0 

680    FOR   1=1   TO   13 
690    C1+C[I3 
700    C2  =  C[I  +  13 
710   Y1=ACI  +  531 
720   Y2=A[I  +  54"3 
730    GOSUB    1160 
740    D1=D1+N1 
750    P3=P3+M1 
760   NEXT    I 
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770  FOR  1=40  TO  52 

780  ci=ccn 

790  C2=C[I+13 

800  Y1=A[I+531 

810  Y2=ACI+54D 

820  GOSUB  1160 

830  D1=D1+N1 

840  P3=P3+M1 

850  NEXT  I 

860  FOR  1=14  TO  39 

870  IF  1=21  THEN  970 

880  C1=-C[I+11 

890  C2=C[I] 

900  Y1=A[I+541 

910  Y2=A[I+533 

920  GOSUB  1160 

930  D2*D2+N1 

940  P4=P4+M1 

950  NEXT  I 

960  GOTO  1200 

970  Cl=-C9 

980  C2=C  21 

990  Y1=Y9 

1000  Y2=A  74 

1010  GOSUB  1160 

1020  D2=D2+N1 

1030  P4=P4+M1 

1040  C1  =  -CC22] 

1050  C2=-C9 

1060  Y1=AL75] 

1070  Y2=Y9 

1080  GOSUB  1160 

1090  D2=D2+N1 

1100  P4=P4+M1 

1110  GOTO  950 

1120  M1=C1*(X2-X1)*(1-0.5*(X2+X1) 

1130  Ml=Ml+(C2-Cl)*(X2-Xl)*(0.5-(X2/3)-(Xl/6)) 

1140  N1=0.5*(C2+C1)*(X2-X1) 

1150  RETURN 

1160  N1=0.5*(C2+C1)*(Y2-Y1) 

1170  M1=0.5*C1*(Y2-Y1)*(Y2+Y1) 

1180  Ml=Ml+(C2-Cl)*(Y2-Yl)*((Yl/6)+(Y2/3)) 

1190  RETURN 

1200  PRINT 

1205  GOTO  Gl  OF  1210,1500,1600 

1210  PRINT  "RUN  NO.  IS"R1 

1220  PRINT 

1225  FIXED  4 

1230  PRINT  "COEFFICIENT  BUILDUP  IS" 

1250  PRINT  TAB2"CNU="L1 

1260  PRINT  TAB2"CNL="L2 

1270  PRINT  TAB10"TOTAL  CN="(L1+L2) 

1280  PRINT  TAB2"CCF="D1 

1290  PRINT  TAB2"CCR="D2 

1300  PRINT  TAB10"TOTAL  CC="(D1+D2) 
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1310  PRINT  TAB2MCMU=MP1 

1320  PRINT  TAB2"CML="P2  » 

1330  PRINT  TAB2"CMF="P3 

1340  PRINT  TAB2MCMR=MP4 

1350  PRINT  TAB10MTOTAL  CM=" (P1+P2+P3+P4) "AT  T.E." 

1360  M1=(P1+P2+P3+P4) 

1370  N1=(L1+L2) 

1380  C1=(D1+D2) 

1390  M5=M1-(0.75*N1) 

1400  DISP  "INPUT  ANGLE  OF  ATTACK  (DEGS.)"; 

1410  INPUT  A2 

1415  A2=0.01745*A2 

1416  A9=A2 

1420  N5=N1*C0S(A2)-C1*SIN(A2) 

1430  C5=C1*C0S(A2)+N1*SIN(A2) 

1440  PRINT 

1445  C2=A2*57.296 

1450  PRINT  TAB  2"CL (MEAN) ="N5 ;TAB21"CD(MEAN) ="C5 ; 

1460  PRINT  TAB2"ALPHA  ="C2"DEGS" 

1470  STANDARD 

1480  GOTO  G3  OF  90,240,250 

1490  END 

1500  N3=L1+L2 

1502  C3=D1+D2 

1504  M3+P1+P2+P3+P4 

1505  FIXED  4 

1506  PRINT  "CN(RMS)*SIN(PHI)="N3 
1508  PRINT  "CC(RMS)*SIN(PHI)="C3 

1510  PRINT  "CM(RMS)*SIN(PHI)="M3 

1511  PRINT 
1595  GOTO  250 
1600  N4=L1+L2 
1602  C4=D1+D2 

1604  M4=P1+P2+P3+P4 

1605  FIXED  4 

1606  PRINT  "CN(RMS)*COS(PHI)="N4 
1608  PRINT  "CC(RMS)*COS(PHI)="C4 

1610  PRINT  "CM(RMS)*COS(PHI)="M4 

1611  PRINT 

1612  Nl=SQR(N3+2+N4+2) 

1613  P1=ATN(N3/N4) 

1614  X=N3 

1615  Y=N4 

1616  F=P1 

1617  GOSUB  1900 

1618  P1=F 

1620  PRINT  "CN(RMS)="N1,"PHI="P1*57.296 
1630  Cl=SQR(C3+2+C4+2) 

1640  P2=ATN(C3/C4) 

1641  X=C3 

1642  Y=C4 

1643  F=P2 

1644  GOSUB  1900 

1645  P2=F 

1650  PRINT  "CC(RMS)="C1,"PHI="P2*57.296 
1660  Ml=SQR(M3+2+M4+2) 
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1670  P3=ATN(M3/M4) 

1671  X=M3 

1672  Y=M4 

1673  F=P3 

1674  GOSUB  1900 

1675  P3=F 

1680  PRINT  "CM(RMS)="M1,"PHI="P3*57.296 
1690  A1=N1*C0S(A9)*SIN(P1)-C1*SIN(A9)*SIN(P2) 
1700  B1=N1*C0S(A9)*C0S(P1)-C1*SIN(A9)*C0S(P2) 
1710  L9=SQR(Al+2+Bl+2) 

1720  L8=ATN(A1/B1) 

1721  X=A1 

1722  Y=B1 

1723  F=L8 

1724  GOSUB  1900 

1725  L8=F 

1730  A2=N1*SIN(A9)*SIN(P1)+C1*C0S(A9)+SIN(P2) 
1740  B2=N1*C0S(P1)*SIN(A9)+C1*C0S(A9)*C0S(P2) 
1750  D9=SQR(A2+2+B2+2) 

1760  D8=ATN(A2/B2) 

1761  X=A2 

1762  Y=B2 

1763  F=D8 

1764  GOSUB  1900 

1765  D8=F 

1770  A3=M1*SIN(P3)-0.75*N1*SIN(P1) 
1780  B3=M1*COS(P3)-0.7  5*N1*COS(P1) 
1790  M9=SQR(A3+2+B3+2) 

1800  M8=ATN(A3/B3) 

1801  X=A3 

1802  Y=B3 

1803  F=M8 

1804  GOSUB  1900 

1805  M8=F 
1810  PRINT 

1820  PRINT  TAB 10 "AERODYNAMIC  COEFFICIENT  SUMMARY" 

1830  PRINT 

1840  PRINT  "CL="N5"  +"L9"SIN(W1*T+"L8*57 . 296") " 

1850  PRINT 

1860  PRINT  "CD="C5"+"D9"SIN(W1*T+"D8*57.296")" 

1870  PRINT 

1880  PRINT  "CM(C/4)="M5"+"M9"SIN(W1*T+"M8*57.296")" 

1890  GOTO  90 

1900  IF  X<0  AND  Y<0  THEN  1920 

1910  IF  X>0  AND  Y<0  THEN  1920 

1915  GOTO  1930 

1920  F=F+PI 

19  30  RETURN 

1940  STOP 
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STEADY  DATA  OUTPUT 


RUN  NO.  IS  123 

COEFFICIENT  BUILDUP  IS 
CNU=  0.3480 
CNL=-0.4281 

TOTAL  CN=-0.0  80  2 
CCF=-0.0339 
CCR=  0.0721 

TOTAL  CC=  0.0382 
CMU=  0.1092 
CML»-0.2840 
CMF*  0.0017 
CMR=  0.0038 

TOTAL  CM=-0.1693   AT  T.E. 


CL(MEAN)=-0.0765   CD(MEAN)=  0.0451    (CM) C/4 (MEAN) =-0 . 1092 
ALPHA  =-4.9991   DEGS 


UNSTEADY  DATA  OUTPUT 


CN (RMS) *SIN (PHI) =-1.538  5 
CC (RMS) *S IN (PHI) =-0.0 31 2 
CM(RMS)*SIN (PHI) =-0.6  336 

CN (RMS) *COS (PHI) =-1.4968 
CC(RMS)*COS(PHI)=-0.004  2 
CM(RMS)*COS (PHI) =-0.8581 

CN(RMS)=  2.1465  PHI=225.7868 

CC(RMS)=  0.0315  PHI=262.2760 

CM(RMS)=  1.0667  PHI=216.4451 


AERODYNAMIC  COEFFICIENT  SUMMARY 

CL=-0.0765  +  2.1405  SIN(W1*T+  225.8305  ) 
CD=  0.0451  +  0.1629  SIN(W1*T+  39.2111  ) 
CM(C/4)=-0.1092   +  0.5836   SIN(W1*T+  63.0435    ) 
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RPPEND1X  B 
PRESSURE  COEFFICIENT  DATA  M1TH  5TEADY  <^ 


RUN  ND.  H2301 

ANGLE  DF  ATTACK=-5  DEB 

Oca0 .  00 
U=100  FT/5EC 

TAP 

PRE55URE 
*     COEFFICIENT 

PRESSURE 
TAP  *     COEFFICIENT 

I 

0.230 

2B 

-0.0B0 

2 

0.537 

29 

-0.115 

3 

0.92B 

30 

-0.0EE 

H 

0.BB7 

31 

-0. 12B 

5 

0.B1E 

32 

-0.107 

E 

0.70E 

33 

-0.133 

7 

0.HH3 

3H 

-0.171 

B 

a.  ire 

3E 

-0.1BE 

3 

0.00E 

3E 

-0.201 

10 

-0.197 

37 

-0.2EI 

11 

-0.259 

3B 

-0 . 2B2 

12 

-0.353 

39 

-0 . 2BH 

13 

-0.H9H 

H0 

-0.339 

1H 

-0.BI I 

HI 

-0.3HB 

IE 

-0.E55 

H2 

-0.391 

IE 

-0.E0H 

H3 

-0.503 

17 

-0.592 

HH 

-0.503 

IB 

-0.591 

HE 

-0.5E9 

19 

-0.HEB 

HE 

-0.E72 

20 

-0.HBE 

H7 

-0.B7E 

21 

-0.390 

HB 

-1.325 

22 

-0.2HH 

H9 

-1 .399 

23 

-0 . 230 

E0 

-0.922 

2H 

-0.075 

51 

-1.H02 

2E 

-0.1HH 

52 

•-0.9H0 

2E 

-0.1EH 

53 

-0.195 

27 

-0.1  IB 
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PRESSURE  COEFFICIENT  PATH  WITH  STEADY  Cy. 


RUN  ND.  H2302 

HNELE  DF  RTTRCK=-5  DEB 


TRP  * 


PRESSURE 
COEFFICIENT 


TRP  * 


=*0.0B55 
100  FT/5EC 

PRESSURE 
COEFFICIENT 


I 

1.000 

2 

0.37H 

3 

0.E39 

H 

0.H7I 

5 

0.233 

E 

0.033 

7 

-0.1BE 

B 

-0.HHB 

9 

-0.E9B 

10 

-0. BIB 

11 

-0.95E 

12 

-1.07H 

13 

-1.205 

1H 

-1 .27H 

IE 

-1.3H3 

IS 

-1.355 

17 

-1.372 

IB 

-1.317 

IS 

- 1 . 235 

20 

-1.190 

21 

-3.371 

22 

-5.552 

23 

-2.B1H 

2H 

-2.HEB 

25 

-1.15H 

2E 

-2.7IB 

27 

-2.9HB 

2B 

0.11E 

29 

-1.171 

30 

-1.072 

31 

-0 . 50H 

32 

-0.379 

33 

-0 . 253 

3H 

0.013 

35 

0.05E 

3E 

0.099 

37 

0.0E1 

3B 

0.023 

39 

-0.105 

H0 

-0. 100 

HI 

-0.09E 

H2 

-0.172 

H3 

-0 . 233 

HH 

-0.221 

H5 

-0.259 

HE 

-0.215 

H7 

-0.279 

HB 

-0.3E3 

H9 

-0.177 

50 

-0.020 

51 

0.!  19 

52 

0.H53 

53 

0.B11 

53 


PRE55URE  COEFFICIENT  ERTR  N1TH  5TERDY  Cg 


RUN  NO.  H2303 

RNELE  DF  RTTRCK=-S  DEE 


TRP  # 


PRE55URE 
COEFFICIENT 


C,u=0.0EH0 
U=100  FT/5EC 

PRE55URE 
TRP  #     COEFFICIENT 


I 

t.000 

2 

0.9B7 

3 

0.7B2 

H 

0.EHI 

E 

0.323 

E 

0.  ISE 

7 

-0.0EI 

B 

-0.3B7 

3 

-0.HBH 

10 

-0.E9S 

11 

-0.B15 

12 

-0.931 

13 

-1.0BE 

14 

-1.150 

IE 

-I.23E 

IE 

-1.23E 

17 

-I. IBS 

IB 

-1.210 

19 

-1.12E 

20 

-1.0EI 

2! 

-2.EB3 

22 

-H.llH 

23 

-I.BB0 

2H 

-0.2HB 

2E 

-0.71B 

25 

-2.0HE 

27 

-1.1SI 

2B 

0.339 

29 

0.203 

30 

-0 . B9E 

31 

-0.279 

32 

-0.22E 

33 

-0.151 

3H 

0.00E 

3E 

0.0IH 

35 

0.023 

37 

-0.023 

3B 

-0 . 03 1 

39 

-0.097 

H0 

-0. I  IB 

HI 

-0.IH0 

H2 

-0.22B 

H3 

-0 . 25B 

HH 

-0.311 

HE 

-0.295 

HE 

-0.313 

H7 

-0.3B7 

HB 

-0.EBH 

H9 

-0.EB7 

E0 

-0.231 

SI 

0.037 

E2 

0.117 

53 

0.709 

54 


FRE55URE  COEFFICIENT  DRTR  N 1 TH  5TERDY  Cu. 


RUN  ND.  4230H 

RNELE  DF  RTTRCK=-5  DEB 


TRP  * 


PRE55URE 
COEFFICIENT 


0**0. 044B 

U=100  FT/5EC 

PRE55URE 
TRP  #     COEFFICIENT 


I 

0.B72 

2 

I.00E 

3 

0.303 

4 

0.EE5 

5 

0.E35 

E 

0.312 

7 

0.030 

B 

-0.197 

9 

-0.3E2 

10 

-0.513 

I  I 

-0.EI9 

12 

-0.775 

13 

-0.92B 

IH 

-1.010 

IE 

-1.092 

IE 

-1.02E 

17 

-1.05H 

IB 

-I. BIB 

13 

-0.BE7 

20 

-0.903 

21 

-1.9H9 

22 

-2.994 

23 

-1.034 

24 

-0.0H5 

25 

-0.229 

25 

-0.7B2 

27 

-0.00B 

2B 

0.237 

29 

0.192 

30 

-0.23B 

31 

-0.0B7 

32 

-0. 13E 

33 

-0.IEH 

34 

-0.0IB 

35 

-0 . 03E 

3E 

-0.05H 

37 

-0.07E 

3B 

-0 . 1 05 

39 

-0.  119 

H0 

-0.143 

HI 

-0. IE7 

H2 

-0.29H 

43 

-0.3E2 

44 

-0.3E7 

45 

-0.3E7 

HE 

-0.3E7 

H7 

-0.534 

HB 

-0.777 

H9 

-0.E13 

50 

-0 . 297 

51 

-0.379 

52 

0.090 

53 

0.579 
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PRESSURE  COEFFICIENT  DRTR  NITH  STEADY  Cu 


RUN  ND .  HH305 

RNELE  DF  RTTRCK=-S  DEE 

O*=0.0272 
U=100  FT/SEC 

TRP 

PRESSURE 
*     COEFFICIENT 

TRP  # 

PRE55URE 
COEFFICIENT 

1 

0.723 

2B 

0.  117 

2 

0.3B0 

23 

0.  I5H 

3 

0.3B3 

30 

0.050 

H 

0.BIB 

31 

0.0H5 

S 

0.E02 

32 

-0.032 

G 

0.H23 

33 

-0.0B2 

7 

0.2EB 

3H 

-0.0H2 

B 

-0.0B7 

35 

-0.0E5 

3 

-0.251 

3E 

-0.0B7 

10 

-0.4  IS 

37 

-0. I2B 

1  1 

-0.503 

3B 

-0. IBE 

12 

-0.E0H 

33 

-0.22H 

13 

-0.75E 

40 

-0 . 22H 

IH 

-0 . B37 

HI 

-0.22H 

15 

-0 . 3 1 B 

42 

-0.230 

IE 

-0.B33 

H3 

-0.33E 

17 

-0 . BEE 

44 

-0 . H32 

IB 

~0.BE! 

HS 

-0.H5I 

13 

-0.E31 

HE 

-0.44B 

20 

-0.B01 

H7 

-0.E03 

21 

-1.370 

HB 

-0.BBB 

22 

- 1 . 3H0 

H3 

-0 . BH2 

23 

-0.53E 

50 

-0 . 743 

2H 

0.057 

51 

-0 . 533 

HS 

-0.0HE 

52 

-2 . 2E2 

HE 

-0.030 

53 

0.  IB3 

27 

0.0EE 

56 


PRESSURE  COEFFICIENT  DRTR  N1TH  5TEHDY  Q. 


RUN  ND.  H230E 

FINELE  OF  RTTRCK=-5  DEE 


TRP  * 


PRESSURE 
COEFFICIENT 


C>»0.0133 
U=100  FT/SEC 

PRESSURE 
TRP  #     COEFFICIENT 


1 

0.522 

2 

0.30a 

3 

0.9EH 

4 

0.B2B 

5 

0.7H3 

E 

0.590 

7 

0.3IB 

B 

0.0HE 

9 

-0.13E 

10 

-0.311 

11 

-0.H07 

12 

-0.HEB 

13 

-0.E0H 

IH 

-0.E97 

IS 

-0.7B9 

IE 

-0 . 7E2 

17 

-0.B50 

IE 

-0.E99 

19 

-0.E32 

20 

-0.EHI 

21 

-0.975 

22 

-1.309 

23 

-0.30E 

2H 

0.00E 

25 

0.0E7 

2E 

0.097 

27 

0.0E7 

2B 

0.0BE 

29 

0.0B5 

30 

0.075 

31 

0.0E3 

32 

0.039 

33 

-0.039 

3H 

-0.0B0 

35 

-0.150 

3E 

-0.220 

37 

-0.279 

3B 

-0.27E 

39 

-0.251 

H0 

-0.291 

HI 

-0.331 

H2 

-0.H07 

H3 

-0.515 

HH 

-0.579 

H5 

-0.5E3 

HE 

-0.E7H 

H7 

-0.772 

HB 

-1 .003 

H9 

-1.123 

50 

-0.H5H 

51 

-1.075 

52 

-2 . B3B 

53 

-0.290 
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RPPENDIX  C 
UNSTEADY  PRE55URE  DRTR 


RNELE  DF  RTTRCK=-5 

DEE          U=100  FT/5EC 

TRP 

MERN  PRE55URE 
#     COEFFICIENT 

UN5TERDY  PRE55URE 
COEFFICIENT 

PHR5E 
RNELE 

1 

0.270 

0.H0H 

1E0 

2 

0.725 

0.1BE 

1E0 

3 

0.30B 

0.073 

230 

H 

0.B12 

0.1EB 

330 

5 

0.750 

0.203 

3H5 

E 

0.E2H 

0.210 

353 

7 

0.3B3 

0.225 

0 

a 

0.1B1 

0.3BH 

17 

a 

-0.0HB 

0.H00 

23 

121 

-0 . 1  BE 

0.5H5 

32 

i  l 

-0.273 

0.513 

3H 

12 

-0.352 

0.EB1 

37 

13 

-0.HB0 

0.730 

HS 

1H 

-0.553 

0.737 

H3 

15 

-0.E2E 

0.7H3 

51 

IE 

-0.E21 

0.7B0 

53 

17 

-0.5B5 

0.B13 

55 

IB 

-0.537 

0.EE1 

E0 

IB 

-0.HBB 

0.500 

EB 

20 

--0.H21 

0.5B2 

75 

2! 

-0.3E7 

0.HE1 

BB 

22 

-0.313 

0.H00 

103 

23 

-0.235 

0.33B 

1  11 

2H 

-0.2B3 

0.2E5 

103 

2S 

-0.275 

0.1B5 

32 

2E 

-0.33B 

0.  IBB 

B3 

27 

-0.H25 

0.  iaa 

E5 

58 


UN5TERDY 

PRESSURE  DRTA 

C CONTINUED) 

BNELE  DF  RTTACK=-5  DEB 

1>100  FT/SEC 

TAP 

MEAN  PRE55URE 
#     COEFFICIENT 

UN5TEADY  PRESSURE 
COEFFICIENT 

PHA5E 
ANGLE 

2B 

-0.H10 

0.  125 

3HB 

29 

-0.2S5 

0.09H 

295 

313 

-0.23B 

0.1E1 

29B 

31 

-0.257 

0.223 

290 

32 

-0.1B1 

0.31H 

290 

33 

-0.173 

0.37E 

2B3 

3H 

-0.17B 

0.H21 

270 

35 

-0.132 

0.H7B 

255 

3B 

-0.205 

0.539 

2H0 

37 

-0.22B 

0.E00 

237 

3B 

-0.2H3 

0.5BB 

230 

39 

-0.2E0 

0.EBB 

231 

H0 

-0.2BH 

0.EBH 

227 

HI 

-0.30B 

0.EB0 

22H 

H2 

-0.3EB 

0.E13 

217 

H3 

-0.50B 

0.S13 

215 

HH 

•-0.55B 

0.E00 

210 

H5 

-0.593 

0.53B 

205 

HE 

-0.EHB 

0.513 

195 

H7 

-0.B53 

0.HE3 

190 

HB 

-1.250 

0.H25 

1B3 

na 

-1 .30B 

0.H05 

175 

52 

-0.9B3 

0.H05 

171 

5! 

-! .453 

0.H50 

152 

52 

-1 .125 

0.50B 

157 

53 

-0.335 

0.H35 

1H2 
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